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Abstract 

Multi-hop relaying is an economically efficient architecture for coverage extension and through- 
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C\J | put enhancement in future wireless networks. OFDM, on the other hand, is a spectrally efficient 
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physical layer modulation technique for broadband transmission. As a natural consequence of 
combining OFDM with multi-hop relaying, the allocation of per-hop subcarrier power and per- 
' hop transmission time is crucial in optimizing the network performance. This paper is concerned 

with the end-to-end information outage in an OFDM based linear relay network. Our goal is to find 
an optimal power and time adaptation policy to minimize the outage probability under a long-term 



total power constraint. We solve the problem in two steps. First, for any given channel realization, 
we derive the minimum short-term power required to meet a target transmission rate. We show that 

it can be obtained through two nested bisection loops. To reduce computational complexity and 

^ ! 

, signalling overhead, we also propose a sub-optimal algorithm. In the second step, we determine a 

m 



power threshold to control the transmission on-off so that the long-term total power constraint is 
satisfied. Numerical examples are provided to illustrate the performance of the proposed power and 



<0 time adaptation schemes with respect to other resource adaptation schemes. 
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I. Introduction 

Relay networks in the form of point-to-multipoint based tree-type or multipoint-to-multipoint 
mesh-type architectures are a promising network topology in future wireless systems. The 
basic concept of relaying is to allow a source node to communicate with a destination node 
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under the help of a single or multiple relay nodes. It has been shown that relaying can bring 
a wireless network various benefits including coverage extension, throughput and system 
capacity enhancement. Recently, multi-hop relaying has been widely adopted in wireless 
networks such as next generation cellular networks, broadband wireless metropolitan area 
networks and wireless local area networks. On the other hand, orthogonal frequency division 
multiplexing (OFDM) is an efficient physical layer modulation technique for broadband 
wireless transmission. It divides the broadband wireless channel into a set of orthogonal 
narrowband subcarriers and hence eliminates the inter-symbol interference. OFDM is one 
of the dominating transmission techniques in many wireless systems, e.g., IEEE 802.16 
(WiMax), EV-DO Revision C and the Long-Term-Evolution (LTE) of UMTS. The combina- 
tion of OFDM and multi-hop relaying has received a lot of attention recently. For example, 
this OFDM-based relay architecture has been proposed by the current wireless standard IEEE 
802. 16] [1]. The complexity of relay station is expected to be much less than the one of legacy 
IEEE 802.16 base stations, thereby reducing infrastructure deployment cost and improving 
the economic viability of IEEE 802.16 systems [2]. 

In this work, we are interested in an OFDM-based linear relay network. The so-called 
linear relay network consists of one-dimensional chain of nodes, including a source node, 
a destination node and several intermediate relay nodes. It can be viewed as an important 
special case of relay networks where only a single route is active. As a natural consequence 
of multi-hop relaying and OFDM transmission, the allocation of per-hop subcarrier power 
and per-hop transmission time is crucial in optimizing the end-to-end network performance. 

Previous work on resource allocation for relay networks is found in [3]-[8]. Yao et al. in 
[3] considered a classic three-node network (a source node, a destination node, and a relay 
node) and compared the energy required for transmitting one information bit in different 
relay protocols. Authors in [4] and [6] studied efficient scheduling and routing schemes in 
one-dimensional multi-hop wireless networks. It is assumed in all these works [3], [4], [6] 
that the point-to-point links are frequency-flat fading channels and the system has a fixed 
short-term power constraint. In [5], Oyman et al. summarized the end-to-end capacity results 
of a multi-hop relay network under fixed-rate and rate-adaptive relaying strategies, and further 
illustrated the merits of multi-hop relaying in cellular mesh networks. Authors in [9] studied 
the per-hop transmission time and subcarrier power allocation problem in the OFDM based 
linear multi-hop relay network to maximize the end-to-end average transmission rate under a 
long-term total power constraint. However, the end-to-end average rate can only be obtained 



at the expense of large delay. 

For many real-time services, one has to maintain a target transmission rate and avoid 
service outage in most of fading condition through adaptive resource allocation. An outage is 
an event that the actual transmission rate is below a prescribed transmission rate ( [10] and 
[11]). Outage probability can be viewed as the fraction of time that a codeword is decoded 
wrongly. For any finite average power constraint, transmission outage may be inevitable over 
fading channels. However, one can minimize the outage probability through adaptive power 
control [10]. In a relay network where no data is allowed to accumulate at any of relay nodes, 
an end-to-end outage is the event that there exists a hop over which the transmission rate is 
lower than the target rate. 

The goal of this paper is to investigate the optimal per-hop power and time control to 
minimize the end-to-end outage probability in an OFDM linear relay network under an 
average transmission power constraint. At first, we derive the minimum short-term power 
required to meet a target transmission rate for any given channel realization. The resulting 
power and time allocation can be obtained through a Two-nested Binary Search (TBS) which 
is conducted in a central controller with the knowledge of channel state information (CSI) 
on all subcarriers and over all hops. Such algorithm gives a theoretical performance limit 
of linear relay networks, but is computationally intense. Moreover, it requires significant 
signalling overhead and channel feedback between network nodes and the central controller. 
For this reason, an Iterative Algorithm of Sub-optimal power and time allocation (IAS) 
is proposed. The required information for signalling exchange only involves the number 
of active subcarriers on each hop and the geometric mean and harmonic mean of channel 
gains averaged over the active subcarriers. This sub-optimal allocation algorithm suggests 
prolonging the transmission time for the hop with low geometric mean of channel gains 
while lowering the transmission power for the hop with low harmonic mean. After obtaining 
the minimum power required to support the target rate for a given channel realization, we then 
compare it with a threshold. The transmission will be cut off if the required minimum total 
power exceeds the threshold. The threshold takes the value so that the long-term total power 
constraint is satisfied. Numerical results show that a significant power saving can be achieved 
by the proposed optimal power and time allocation compared with the uniform power and 
time allocation under the same end-to-end outage probability. In addition, the proposed sub- 
optimal power and time allocation serves as a good approximation to the optimal solution 
when the target rate is sufficiently high. The optimal number of hops in the sense of requiring 



minimum power at different target rates is also shown numerically. 

The remainder of this paper is organized as follows. In Section [II], the system model 
and problem formulation are presented. The optimal and sub-optimal resource allocation 
algorithms to minimize the end-to-end outage probability under an average total power 
constraint are proposed in Section |nij Numerical results are given in Section [IV] Finally 
we conclude this paper in Section |V] 

II. System model, end-to-end rate and outage probability 

A. System Model 

Consider a wireless linear relay network shown in Fig. [H The source node R commu- 
nicates with destination R N by routing its data through N — 1 intermediate relay nodes 
Rn (n — 1, . . . , N — 1). The hop between node R n -i and R n is indexed by n, and the set 
of hops is denoted by J\f. We focus on time-division based relaying. The transmission time 
is divided into frames of multiple time slots. Within every time frame, the transmission over 
each hop takes place at the assigned time slots. In general, frequency reuse can be applied so 
that more than one hop can be transmitting at a same time slot. However, due to interference 
issue, it will increase decoding complexity as well as decoding delay [4]. Thus, in this work 
we do not pursue the frequency reuse. In each time frame, the message from the source is 
sequentially relayed at each hop using decode-and-forward protocol [12]. Each relay decodes 
the message forwarded by the previous node, re-encodes it, and then transmits it to the next 
receiver. The channel for each hop is assumed to be a block fading Gaussian channel, and the 
channel coefficients remain constant during the entire time frame but change randomly from 
one frame to another. Over each hop, OFDM with K subcarriers is used as the physical layer 
modulation scheme. We denote the set of subcarriers by /C. The channel gain on subcarrier 
k over hop n in a time frame is denoted as and it is independent for different n. 

B. End-to-end rate and outage probability 

Suppose each time frame contains T OFDM symbols and hop n is scheduled to transmit 
over T n OFDM symbols with T n satisfying Yln&M Tn = T. Then we define the time-sharing 
fraction as p n = 2*. It is assumed that T is large enough so that p n can take an arbitrary 
value between and 1. Let pk : n denote the transmission power on subcarrier k over hop n. 



It is subject to a long-term total power constraint P, given by 



E 



< P- (1) 



_neN k&K. 

Then the instantaneous transmission rate in Nat/OFDM symbol in a time frame achieved 
over hop n can be written as 



rn = p„$> 1 + , Vn E A/", (2) 

where iVo is the noise power, and Y is the signal-to-noise ratio (SNR) gap related to a given 
bit-error-rate (BER) constraint [13]. For notation brevity, in the remaining part of this paper, 
we redefine g^n as g kn := gk,n/(TN ). Under the assumption that no data is allowed to 
accumulate at any relay nodes (also called "information-continuous relaying" in [5]), the 
total number of bits received at the destination node at the end of time frame, B, is the 
minimum of the number of bits transmitted over each hop, B n , where B n = r n T n . Thus, 
the end-to-end transmission rate r can be defined as r = min neA rr n . In the following we 
introduce the end-to-end rates under different resource adaptation policies. 

Uniform power and time allocation (UPT): When each transmitting node has no CSI, or 
does not exploit CSI due to high signalling overhead, the transmission scheme is independent 
of the CSI and both the time and power are uniformly allocated. Hence, the end-to-end rate 
can be achieved as 

where g = {gk, n , k E K,n E A/"}. As can be seen, the end-to-end rate is limited by the hop 
with the worst channel condition. 

Fixed power and adaptive time allocation (FPAT): When the transmitters have CSI to 
some extend (not necessarily global CSI), each node can perform rate adaptation to avoid the 
situation that the ill-conditioned hop become the bottleneck of the whole link. We assume 
that the transmission power on each subcarrer over each hop keeps unchanged, and rate- 
adaptation is performed by adjusting time-sharing fraction such that r n = rj,Vi ^ n E N. 
In this scenario, the maximum end-to-end transmission rate is given by [14] 

r(g '"' F) = fe um'^ffl ) ■ <4) 

This rate is achieved by assigning the time sharing fractions p = {p. h i E A/"} to be 
f , IIn*E fce K [^(l + g k , n P/K)] + . 

EnGAT Uj^n EfcgK (1 + 9kjP/K)\ 







By comparing (HJ) with ©, it is found that the end-to-end rate is increased by adaptive time 
allocation. This is because the harmonic mean of a set of nonnegative values is always greater 
than or equal to the minimum value. To implement FPTA, the central controller only needs 
each hop to feedback the value of £] fe6;c ln(l + 9k, n P/K) instead of collecting the global 



Adaptive power and fixed time allocation (APFT): In this case the time-sharing fractions 
are assumed to be fixed and equal to each other, but the transmission power can be adjusted 
adaptively. Hence, the conditional end-to-end rate for a given power value set p = {pk, n , k G 
JC,n E A/"} is expressed as 



The set of power values p depends on the global CSI g and the target end-to-end rate. 

Adaptive power and time allocation (APT): We shall now focus on the scenario of interest, 
where both transmission power over each subcarrier and each hop and time over each hop 
are allowed to be dynamically allocated. We assume that at the start of each time frame, the 
global CSI is perfectly known at a central controller, which could be embedded in the source 
node. The instantaneous end-to-end rate for given power values p and time-sharing fractions 
p is expressed as 



L fce/c J 
Let p satisfy the time constraint YlmeN P n = an ^ P satisfy the long-term power constraint 
([U). The end-to-end information outage probability evaluated at target rate R for APT can 
be expressed as 



Our goal is to minimize P^p t T (R, P) with respect to the power and time adaption {p(g), p(g)}- 
Namely, 



CSI g. 




(6) 




(7) 



P$ T (R,P) = P(r(g,p,p)<R). 



(8) 



PI: min 

{Pn ,Pk,n} 
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The next section is dedicated to solving the problem PI. As it will be clear later, APFT 
can be viewed as a special case of APT by fixing p n = 1/N, Vrt and hence the minimization 
of its outage probability can be solved similarly. 



The minimum outage probability problem PI defined in the previous section can be 
generally solved in two steps as proposed in [10]. First, for each global channel state g, 
the short-term minimum total power p m i n (g) required to guarantee the target end-to-end 
transmission rate R is to be determined. The second step then determines a threshold to 
control the transmission on-off subject to a long-term power constraint. 

A. Short-Term Power Minimization 

In this subsection, we shall find the optimal time sharing fraction p* (Vn G A/") and the 
optimal power allocation p* k n (Vrt G M, k G K) to minimize the short-term total power 
needed to achieve a target end-to-end transmission rate R. Then a sub-optimal algorithm 
with reduced complexity is developed. The sub-optimal one has a closed-form expression 
from which a few attractive properties regarding time and power allocation can be observed. 
Comparison on average powers and computational complexity between the optimal and sub- 
optimal algorithms is also given. 

1) Optimal power and time allocation: The optimal power and time allocation problem 
to minimize short-term total power can be formulated as 



Unfortunately, the function r(g, p, p) defined in © is not concave in p and p. As a result, 
the problem P2 is not convex. To make the problem P2 more tractable, we introduce a 
new variable s fe)Tl defined as Sk, n '■= PnPk,n- This new variable can be viewed as the actual 
amount of energy consumed by hop n on subcarrier A; in a time frame interval. In addition, it 
follows from © that constraint dT2]) can be rewritten as N sub-constraints. By doing these, 
problem P2 is transformed into a new problem with optimization variables p n (Vn G Af) and 



III. Adaptive power and time allocation (APT) 
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s.t. r(g, p, p) > R 



(12) 




Sk, n (Vn e Af, keJC): 



P3: min 

{pn, Sk,n} 



s.t. 
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Since its Hessian matrix is negative semidefinite, the function p n ln(l+g kj1l Sk , n / Pn) is concave 
in p n and st, n - Therefore, problem P3 is a convex optimization problem and there exists 
a unique optimal solution. To observe the structure of the optimal solution, we write the 
Lagrangian of Problem P3 as follows: 

J({Pn},{ S k,n},{ X n},P) = ^2^2 s k,n + P I ^ Pn ~ 1 I + 

neAf keK \n£j\f ) 



i?-p n Vln ( 1 + 9J^A 



(16) 



where A„ > (n G A/") and /3 > are the Lagrange multipliers for the constraints (TT4l) 
and (fT5l) . respectively. If {p*} and {s£ n } are the optimal solution of P3, they should satisfy 
the Karush-Kuhn-Tucker (KKT) conditions [15], which are necessary and sufficient for the 
optimality. The KKT conditions are listed as follows: 

dJ{. . .) ^ 

f *tn = 



ifst n >0 



, Vn G Af, k G K 



(17) 
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dJ{...) 



dp n 



>o if P ; = o 

= if < p* n < 1 , Vn G TV 

<o if P ; = i 
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Pi 



R 



0, Vn G Af. 



(18) 



(19) 



It can be obtained from the KKT condition (flTl) that the optimal power distribution {p* k n } 
has a water-filling structure, and is given by 



Pk,- 



>k,n 



A, 



VfcG/C, n £ Af, 



(20) 



Pn V #fc,n, 

where (a;) + = max(0, x), and A n can be regarded as the water level on hop n. Different hops 
may have different water levels. For each hop, more power is allocated to the subcarrier with 
higher channel gain and vice versa. 



Let K n denote the set of subcarriers over hop n that are assigned with non-zero power, 
i.e. satisfying g k n > 1/A n , VA; G K, n , and let k n be the size of the set. The subcarriers in 
the set are said to be active subcarriers. Note that each water level value A n cannot be zero. 
Otherwise p* k n = 0, Vfc, n and, as a result, the constraint (TT4l) cannot be satisfied. Hence, we 
obtain the closed-form expression for p* by substituting (l20l) into the KKT condition ( TT9T ): 



Pn = US, K) = ^ : ^— T-j-p , Vn G JV. (21) 

From (12TT) . it can be shown that p* is monotonically decreasing in A n (note that k n also 
depends on A n ). 

In the following, we derive the relation between \ n and (3. Taking the derivative of 
Lagrangian of P3 in (fT6l) with respect to p n , we have 



dJ{...) _ 

— A T 



9k,n$k,n \ \ ^ 9k,n^k,n 



Pn J keK . Pn + 9k,nS k ,: 



- P- (22) 



dp n 

Suppose that there exists an n G M such that p* = or 1, then the condition (fT9l) would be 
violated. Thus, we have < p* < 1 , Vn G N . Substituting ([20]) into d22j) and using (fT8l) . 
we express /5 as a function of A ra given by 



/5 - /n(g, A n ) = A n ( ^ lng- fc , n + fc n In A n - k n | + ^ — , Wn G A/". 



(23) 



It is seen from (|23T) that finding the optimal water levels {A n } at a given (3 are independent 
problems. It can be proven that f n (g, A n ) is a monotonically increasing function of A n in 



the region min^ ( ) , +oo by evaluating its derivative with respect to A n . Hence, the 

V 9k, n J 

inverse function, A n = /~ 1 (g, (3), exists and is an increasing function of f3. Therefore, the 
exact value of A^ for a given (3 can be obtained numerically using binary search when the 
upper bound is known. 

Substituting A n = /^(g,/?) into dlB, we can express p* n as p* n = h^f^fe, f3)). Since p* 
is decreasing in A n and A n is increasing in (3, we have that p* is decreasing in (3. Therefore, 
the optimal (3 can also be obtained via binary search from the constraint (TT5T) . Hence, the 
optimal solution {p* , s* k n } of P3 and the resulting p min can be obtained through two nested 
binary searching loops. The outer loop varies the Lagrange multiplier (3 to meet the time 
constraint. The inner loop searches the water level for each hop at a given value of (3 to 
satisfy (l23l) . The algorithm is outlined as follows. 



LU 



Two-nested Binary Search for minimum short-term power (TBS) 
Binary search for /3 

1) Find the upper bound and lower bound of (3 

a) For all n G TV", let A„ = A„ = vain k {l / g k , n } 

b) Compute p n = h n (g,X n ) using (l2TT i 

c) If p n > 1/N, update A n = 2A n and repeat Step l)-b) and c), else, go to Step l)-d) 

d) Set /3 min = max„ ejV /„(g, AJ, and /3 max = max„ eAf / n (g, A„) 

2) Set high = /3 max , low = /3 min 

3) Let center = (low + high)/2 and binary search for A n (Vn e Af) at /3 = center 

a) Find the upper bound and lower bound of A ra , A™ ax and A™ m , respectively 

i) Let A£ in = A" ax = min fe {l/ 5fe , n } 

ii) Compute (3' = ./r>(g,A™ ax ) using <|23]> 

iii) If /?' < /3, update A™ 11 = A„ lax and A^ lax = 2A£ lax and repeat Step 3)-a)-ii) and in) 
else, let high n = A'f ax , low n = A™ in , and go to Step 3)-b) 

b) Set center n = (low n + high n )/2. If / n (g, center n ) > (3, let high n = center n ; otherwise, let 
low n = center n 

c) Repeat Step 3)-b) until high n — low n < e' and let A„ = center n 

4) If J^neAf hn(g, A n ) > 1, let /ow = center; otherwise, let high = center 

5) Repeat Step 3) and Step 4) until high — low < e 

6) Using the found {A„} and (3, obtain p* and p* k based on (fJTJ and d20l >. respectively. 

7) Compute Pmiu = EngJVPnCEfeggPfe.n) 

In Step 1), the boundaries of (3 are determined in order to proceed with the binary search 
in the outer loop. From (1231) . a common Lagrange multiplier (3 is shared by all hops and it is 
a monotonically increasing function of A n in the region of minfcgA: ( — ^— ) , +oo for all n. 
We use A = min fcg ^ ( ) to represent the lower bound of A n , then the lower bound of (3 
is maximum of / n (g, A„) among hops. The same lower bound of A n will also be used in 
Step 3) for the inner loop. The upper bound of j3 is obtained from the fact that there exists at 
least an n* such that p n * > jv. Correspondingly, we find a water level A n = h~ 1 (g, ■h) for all 
n, where /i^ 1 (g, •) is the inverse function of h n (g, •). Then for hop n*, we have A n * < A n * 
since h~ 1 (g, •) is a decreasing function. Therefore, due to the monotone of / n (g, •), the upper 
bound of (3 can be obtained from 

P = Ms, V) < Ms, V) < max/ n (g, A n ). 
The algorithm then updates j3 using binary chop until the sum of the corresponding time- 
sharing fraction converges to 1 . The convergence of the outer loop is guaranteed by the fact 



that the actual sum of time-sharing fractions is also monotonically decreasing in (3. 

The aim of the inner loop in Step 3) is to find A n (Vrt G TV) for a given (3. We first 
initialize the upper bound of A n , A™ ax = min^ j^— j and then keep increasing it until the 
corresponding (3' goes beyond the given (3. In each iteration, the binary search guesses an 
halfway A n between the new high and low and repeats it until the actual f3' approach the 
given (3. The iteration converges because of the monotone of (3 in A n . 

The outer loop involves log 2 ^ / 3max ~/ gmin j iterations where e represents outer loop accuracy. 
The inner loop has N binary searches, and each involves log 2 ( - -, - ) iterations, where e 
is the inter loop accuracy. It is observed from (ED and ([23]) that (3 max = O (NRe NR/K ) and 
A™ ax = O (e NR / K ^ when the target rate is so high that all subcarriers are active. Therefore, 
the average computational complexity of this algorithm is upper bounded by the magnitude 
of n ^ r2 ln(^-) ln(p-) in the asymptotic sense with a high target rate. 

2) Sub-optimal time and power allocation: In the optimal time and power allocation, 
it is infeasible to obtain an closed-form expression for the solution. In this part, we will 
observe that when the target rate is sufficiently large, the optimal transmission time can be 
approximated by an explicit function of geometric mean of channel gains averaged on active 
subcarriers and the number of active subcarriers. In addition, the product of water level and 
the number of active subcarriers for each hop tends to be the same. In the following, we shall 
investigate this sub-optimal solution and show that it has a low computational complexity 
and little signalling exchange. 

Let {p n , n G J\f} be any given time allocation that satisfies J2neAf Pn = 1 and < p n < 1. 
The optimal power distribution at the given {p n , n G A/"} is expressed by (l20l) . Substituting 
(|2Q|) into © and letting p n = Ylk&K.Pk,n, the close-form expression of A„ can be obtained 



where k n is the size of the set /C n = {k \ g kn > 1/A n }. Moreover, substituting (l24l) back 
into (T20l) . we have 



as [16] 




(24) 



Pn(pn) 




(25) 



where, for notation brevity, we define 




(26) 



and 

k£lCn 9k,n 9n 

In (1261) and (l27l) . g n and represent the geometric mean and harmonic mean of gk, n over 
the active subcarriers at hop n, respectively. 

For the moment, let us assume that &; n 's are fixed and then both a n and b n are constants. 
Then, the problem of minimizing total power for supporting the target end-to-end transmission 
rate can be reformulated as P4 only with optimization variables {p n , n G M} 

P4: Pmin(g) = min V] p„p* = min V" p n (e^^~ a " - 6„ ) (28) 
s.t.^ Pn = l. (29) 

n£Af 

Problem P4 can be also solved using Lagrange multiplier method since it is convex. The 
Lagrangian of this problem is given by 

L( P , 1S) = J2 (Pn^^~ an ~ &„) + V I 1 - Pn ) , 



where the lagrange multiplier v satisfies constraint (1291) . Applying KKT condition, the optimal 
time-sharing fraction p n should satisfy 

dL(p.u) _r_j_ „ R 1 _rj_ „ 

^ ' = e fc « pn ~ an - - e fe n pn ~ an - v = o. (30) 

Opn fan Pn 

The closed-form solution to (|3Q|) is difficult to obtain. 

It is known that when the target transmission rate is sufficiently small, the power saving 
through time adaptation is insignificant [17]. This result motivates us to focus on the high 
target transmission rate with R 3> K. We consider two particular hops, ri\ and n 2 . Under 
the assumption of a high target transmission rate, the equation (|30l) can be approximated by 

R s _ R R 

^niPni k n2 p n2 

From the above approximation, we can obtain a ratio 

kmpri! 1 k n p ni k ni p ni I k ni p ni \ 

H d — i a n 2 - a„J — In . (31) 



kn,2pn2 R R \kri2Pn2 / 

Without loss of generality, we assume a ri2 > a ni , then we have (k ni p ni )/(k n2 p n2 ) > 1 from 
(|31) . Thus, (HB leads to 

l<|^<l + ^^K 2 -a ni ). (32) 

^ri2Pri2 



Using the inequality ln(x) < x — 1 and (13TI) . we obtain a lower bound of (k ni p ni ) / (k m p n2 ) 
after some manipulations, 

fcniPni ^ ^ _|_ k ni p ni (a n2 — a Wl ) (33) 
Since i? 3> K , inequalities (|32|) and (|33T) lead to 

kniPni ^ k ni p ni 
kn2Pn,2 

After manipulation, we have 

1 CLno 1 ^ni w » r 

— i Vn 1 ,n 2 e A/. 



kn2pn2 H k ni p ni R 
Let 

/i, Vn G A/". 



We can obtain the approximated but close-form solution to (l30l) as follows 



k n Rp + k n CL n 



where /i is determined by the time constraint (1291) . and can be obtained through binary search. 
Substituting (|34|) into (|25T) . the corresponding transmission power allocated to hop n is given 
by 

p' n = e* i -b n . (35) 



Furthermore, substituting (1341) into (|24|) yields the sub-optimal water level for hop n as 

K = ~r ■ (36) 

From (|2Q|) . (1351) and (l36l) . we can regard the sub-optimal power allocation algorithm as a 
two-level water filling algorithm. First, the power is poured among all the hop according to 
(1351) using the water level e R/i and the hop with small b n will be given more power. In each 
hop, the power obtained from the previous level is then poured among different subcarriers 
following (|2Q|) . and the water level is equal to e R/i //c n . 

Consider a special case where e Rfl is sufficiently large so that all subcarriers are active, 
i.e., k n = K. It follows immediately from (|34|) that the hops with low geometric mean of 
channel gains over the subcarriers should be assigned with longer transmission time. Also, it 
follows from (|35T) one should lower the transmission power for the hops with low harmonic 
mean of channel gains. An intuitive understanding of this result is that a high priority is given 
the hop with poor channel condition to take advantage of "Lazy Scheduling" [18] to avoid 



this hop becoming the bottleneck of the whole link. The idea behind "Lazy Scheduling" is 
that energy required to transmit a certain amount of information decrease when prolonging 
transmission time. 

We now relax the assumption made earlier that A; n 's fixed and propose an iterative procedure 
to find the best k n 's for this sub-optimal problem. 

Iterative Algorithm of Sub-optimal power and time allocation (IAS) 

1) Initialization of k n 
Set kn = K, Vn G TV 

2) Binary search for \i for a given {k n , Vn} 

a) Set high = ^ max , low = ^ min 

b) Let center = (low + high)/ 2, and calculate {p' n , Vn} when fi = center according to (l34l 

c ) H ^2 nG j\f p'n > 1> l et l° w ~ center; otherwise, let high = center 

d) Repeat Step 2)-b) and c) until high — low < e" 

3) Find k„ (Vn) in the set {1, . . . , K} for a given p' n to meet the target rate R based on (l37b 

4) Repeat Step 2) and 3) until k n 's are unchanged 

5) Compute p' k n through substituting (l24l into (|20| > 

6) Obtain the required total power p^ in = ^(Efc Pfe,„) 

In Step 2)-a) /i max and /i m i n represent the upper bound and lower bound of fi, respectively. 
The exact value /i min = max ng ^ ^ ^ — ^-J can be obtain from the time constraint < p' n < 
1. Its upper bound // mEPC could be min neA r since if fj, > mm neAf O^jf 21 \ 

p' n < l/N (Vn) from (|34|) . which violates the constraint J2 n Pn = 1- 

The implementation of the above algorithm can be done as follows. At the beginning of 
each time frame, we first assume that the transmission is on for all subcarriers. The central 
controller searches for jx and broadcast it to all relays. The relays and source node then 
compute their own transmission time {p' n , Vn e M} using (|34|) locally. The required power 
allocation for hop n to meet the target rate should satisfy 

^2 m (! + 9k,nPk,n) = ^2 ln (9k,nK) = — , Vn E jV. (37) 

The left side of the above equation (a) can be shown to be a monotonically increasing function 
of A n , and is denoted as z n (X n ). With loss of generality, we assume g l n > g 2>n > . . . > 
9k n (Vn e /C). Each A n maps to a unique k n which satisfies that < \ n < — - — . Thus, 



we have 



R 



^9k n ,n J P n \9k n +l,n. 

Therefore, the desired k n in Step 3) can be obtain through binary search in the set of 
{1,...,K} by comparing z n {l/g kn , n ) with Rj ' p' n . The found k n and the geometric and 
harmonic mean of channel gains on these k n subcarriers are returned to the input of (l34l) 
in the central controller. This procedure repeats until the k n 's, are unchanged. Although the 
convergence of this algorithm cannot be guaranteed theoretically, divergent behaviors were 
never observed in the simulation. In the following, we shall use simulation to examine the 
average number of iterations for the algorithm to converge and the average required short-term 
total transmission power. 

In the simulation, SUI channel model for the fixed broadband wireless access channel 
environments [19] is used and the channel parameters will be detailed in Section |IV] The 
simulation is run for 10 3 time frames to evaluate the average performance. The number of 
subcarriers is set to 16. 

Fig. [2] shows the average iterations in the outer loop over 10 3 independent channel realiza- 
tions required for the search of {k n } to converge. It is shown that the average iteration 
numbers, denoted as M, is decreasing in R and approaches 1 when the target rate is 
sufficiently large. It can be explained by the fact that k n = K when R goes infinity. 

Since the binary search for /x in the inner loop involves log 2 ( /Xmax ~ / Atmin ) iterations and 
finding k n for a given Rj p n involves log 2 (i^) ones, the total number of iterations required 
for the IAS can be express as 



C IAS = M { log 2 



+ iVlog 2 (/m (38) 



Since M is decreasing in R, Cias is also decreasing in R and upper bounded by a linear 
function of N. Fig. [3] compares average total complexities between TBS and IAS for different 
R and different N. 

Fig. |4] compares the average power required to meet the target rate between TBS developed 
in Section Ull-A. 1 1 and its sub-optimal algorithm, IAS. It is shown that IAS serves as a good 
approximation of TBS, especially for a high target rate. 

As we discuss previously, the required controlling signals from the feedback channel are 
only geometric mean, harmonic mean of gk, n and the number of active subcarriers over each 
hop instead of {gk, n ,Vfc G K,,n E M} as in TBS, thus the signalling exchange is greatly 



reduced when the number of subcarriers is large and/or the target rate is high. Furthermore, 
since IAS has low complexity and near-optimal power consumption performance, it is a good 
candidate for a sufficiently high target rate in a real system. 

B. Long-Term Power Threshold Determination 

We have discussed the short-term total transmission power minimization. If the transmis- 
sion is on for every possible channel realization, the long term power constraint may be 
violated. Similar to the single user case [10], the optimal power allocated to all hops for PI 
with a long term power constraint must have the following structure, 

, , I Pmin(g) with probability w(g) 
Kg) = S • (39) 

with probability 1 — w(g) 

Thus, the outage probability is P(r(g, p, p) < R) = E[l — w(g)]. Then solving PI is 
equivalent to finding the optimal weighting function w(g) to the following problem, 

min E[l — u>(g)] 
s.t. < w(g) < 1 

E[ Pmin (g)w(g)} = P. 

According to the result of [10, Lemma 3], the optimal weighting function has the form 

1 for p m in(g) < S* 

w*(g)={ w for p min (g) = s* . (40) 
for p min (g) > s* 

The power threshold s* is given by s* = sup{s : V(s) < P}, and w is given by w = 
where the re gi on n ( s ) and K(s) are defined as K(s) = {g : p min (g) < 
s}, TZ(s) = {g : Pmin(g) < s}, and the corresponding average power over the two sets are: 
V(s) = E ge ft (s) [p min (g)],P(s) = E ge7 ^ (s) [j9 min (g)] The resulting minimum outage probability 
is denoted as 

P£pt = 1 " ProMg G TZ(s*)} - w Pwb{p min (s) = s}. 



From (1391) and (1401) . we see that when the minimum total power for all hops required to 
support the target transmission rate is beyond the threshold s*, transmission is turned off. 
When the required power is less than the threshold, the transmission follows the minimum 
transmission power strategy derived from Section IIII-AI 



The value of s* can be computed a priori if the fading statistics are known. Otherwise, 
the threshold can be estimated using fading samples. During the estimation of the threshold, 
since the channel is assumed to be ergodic, the ensemble average transmission power is equal 
to the time average 

1 * 

E gG TC( s )[Pmin(g)] = Jim - 

t — >oo t — 

i=l 

where pt^ represents the actual transmission power at time frame i. Thus, the threshold is 
always adjusted in the opposite direction of P — | Yll=i P(i) as 



s*(t 



s*(t) 



1 + e [P- 



i=i / 



(41) 



where t is the time frame index. 

Combining the short-term power minimization and long-term power threshold determina- 
tion, the full algorithm for APT is outlined as follows. 



APT 



1) Set t = 1 and s* t) = P 

2) Search for minimum short-term power (developed in Section IIII-AI) 

3) On-off decision 

If Pmin > s%y turn off the transmission and let p( t ) = 0; otherwise, turn on the transmission and let 

P(t) =iW- 

4) Update the threshold s* 



s*(t+l) = s*(t) 1 + e 
5) Let t = t + 1 and return to Step 1). 



i=l 



(42) 



If TBS in Section IIII-A.ll is used in Step 2), we name the optimal APT as APT-opt for 
short. If IAS in Section IIII-A.2I is used, we denote it as APT-sub. 



C. Special case: adaptive power and fixed time allocation 

APFT can be viewed as a special case of APT by fixing p n = 1/N. It can also be solved 
following two steps: short-term power minimization and long-term power threshold determi- 
nation. However, unlike APT, the first step can be performed locally, i.e., each transmitter 



only needs to know the local CSI over the associated hop to solve the problem 

min p kn 

S-t. ^ M 1 + 9k,nPk,n) > R, 

keic 

for all n E M. The solution of the problem is easily obtained as (|2Q|) . where the water level 
is given by ([24]) with p n = 1 /N. 

IV. Numerical results 

In this section, we present some numerical results to illustrate the performance of the 
proposed adaptive power and time allocation for OFDM based linear relay networks. The 
proposed algorithms, APT-opt and APT-sub, are compared with UPT, FPAT and APFT as 
defined in Section HB 

We consider an iV-hop linear wireless network. The acceptable BER is chosen to be 1CT 5 , 
which corresponds to 8.2dB SNR gap. We fix the bandwidth to be 1MHz and the end- 
to-end distance to be 1km. The relays are equally spaced. In all simulations, the channel 
over each hop is modelled by Stanford University Interim (SUI)-3 channel model with a 
central frequency at around 1.9 GHz to simulate the fixed broadband wireless access channel 
environments [19]. The SUI-3 channel is a 3-tap channel. The received signal fading on the 
first tap is characterized by a Ricean distribution with K-factor equal to 1. The fading on 
the other two taps follows a Rayleigh distribution. The root-mean- square (rms) delay spread 
is 0.305/iS. Then the coherence bandwidth is approximately 65KHz. Hence, the number of 
subcarrier K should be greater than 15.2 so that the subcarrier bandwidth is small enough to 
experience the flat fading. Here we choose K = 16. Doppler maximal frequency is set to 0.4 
Hz. Intermediate path loss condition ( [20, Category B]) is chosen as the path loss model, 
which is given by PL = A + a lg J , where A = 20 lg(4.7rdo/^) (A being the wavelength 
in m), a is the path-loss exponent with a = (a — bhb + c/hb). Here hb = 30m is chosen 
as the height of the base station , do = 100m and a, b, c are 4, 0.0065 and 17.1 given in 
[20]. The corresponding a will be used in all simulations except the one in Fig. [6[ In each 
simulation, 10 4 time frames are used to estimate the outage probability. 

Fig. [5] shows the end-to-end outage probabilities versus average total transmission power 
for R = 1,20 and 40 Nat/OFDM symbol using APT-opt when N varies in the set of 
{1,3,5}. From the figure, it is shown that multi-hop transmission can help to save total 



power consumption when the target transmission rate is low (e.g., R = 1) whereas it is better 
to send data directly to the destination if the target transmission rate is high (e.g., R = 40). 
That can be explained by the following fact. As the number of hops increases, the path loss 
attenuation on each hop reduces. But the transmission time spent at each hop also reduces 
since the total frame length is fixed. It is observed from © that the transmission rate is 
linear in transmission time and concave in channel gain. Hence, when the target transmission 
rate increases, the loss due to transmission time reduction cannot be evened out by the gain 
brought by path loss reduction. 

Fig. [6] shows the optimal number of hops to achieve minimum power consumption at 
different target transmission rates. Here, the outage probability is fixed to 1%, and the path 
loss exponent a =2.5 and 4, respectively. It is observed that the optimal number of hops 
is roughly proportional to the inverse of R, and increasing linearly in a. A similar trend is 
shown in [21] where a spacial case, frequency-flat fading channel and a fixed short-term total 
power constraint, is considered. 

Fig. U\ and Fig. [8] compare the end-to-end outage probabilities achieved by different power 
and time adaptation schemes for R = 1 and 20 Nat/OFDM symbols, respectively. A number 
of interesting observations can be made from the two figures. First, by comparing the curves 
of FPAT and UPT it is observed that just adapting per-hop transmission time alone can 
increase the performance considerably. But the decreasing speed of the outage probability 
as the total power increases is not increased much. On the other hand, by comparing the 
curves of APFT and UPT, it is seen that power adaption can bring dramatic improvement on 
the performance. In particular, the slope of the outage probability curves approaches almost 
infinity. This indicates that by turning off the transmission when the channel suffers from 
deep fade can achieve significant power saving. Next, comparing APT-opt with APFT we 
can see that time adaptation on top of power adaptation is still beneficial, but the gain is 
rather limited when the target data rate is small. Finally, it can be seen that the performance 
of APT-sub is even worse than that of APFT when the target rate is low (e.g. R = 1). But 
for large target rate (R = 20), APT-sub becomes superior and is near optimal. 

The above numerical results suggest that multi-hop transmission is favorable at low and 
medium target rates, whereas a direct transmission from source to destination is preferred 
if the target rate is high. Also, power adaptation plays a more important rule than time 
adaptation in minimizing the end-to-end outage probability. In particular, APFT is a good 
choice in practice for low target rates since it has similar performance with APT-opt and yet 



zu 

is much less complex. For the similar reason, APT-sub is recommended at medium target 
rates. 

V. Conclusions 

In this work, we consider adaptive power and time allocations for OFDM based linear 
relay networks for end-to-end outage probability minimization. The problem is solved in two 
steps. First, we derive the minimum short-term total power to meet the target transmission 
rate. Both optimal and sub-optimal algorithms are proposed. In particular, the sub-optimal 
algorithm suggests prolonging the transmission time for the hop with low geometric mean 
of channel gains averaged over subcarriers while lowering the transmission power for the 
hop with low harmonic mean. In the second step, the transmission on-off is determined by 
comparing the required minimum total power with a threshold, which is selected to satisfy the 
long-term total power constraint. Numerical study is carried out to illustrate the performance 
of different resource adaptation schemes: APT-opt, APT-sub, APFT, FPAT and UPT. We find 
that the three schemes with adaptive power control, APT-opt, APT-sub and APFT, provide 
significant power savings at a same end-to-end outage probability over the other two. While 
APFT is a good choice for practical implementation at low target rates, APT-sub becomes 
near optimal at medium target rats. 
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Fig. 2. Average number of iterations in the outer loop required for the search of {k„} 
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Fig. 3. Average total number of iterations using TBS and IAS 
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Fig. 5. End-to-end outage probability vs. average total transmission power under APTA when K — 16 
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Fig. 6. The optimal number of hops vs. target rate under APTA when a = 2.5 and 4 
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Fig. 7. End-to-end outage probability vs. average total transmission power under APT-opt, APT-sub, APFT, FPAT and 
UPT when K = 16 and R = 1 Nat/OFDM symbol 




Fig. 8. End-to-end outage probability vs. average total transmission power under APT-opt, APT-sub, APFT, FPAT and 
UPT when K = 16 and R = 20 Nat/OFDM symbol 



